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Abstract 

We adopt a covariant formalism to derive exact evolution equations for nonlinear perturbations, in 
a universe dominated by two scalar fields. These scalar fields are characterized by non-canonical 
kinetic terms and an arbitrary field space metric, a situation typically encountered in inflationary 
models inspired by string theory. We decompose the nonlinear scalar perturbations into adiabatic 
and entropy modes, generalizing the definition adopted in the linear theory, and we derive the 
corresponding exact evolution equations. We also obtain a nonlinear generalization of the curvature 
perturbation on uniform density hypersurfaces, showing that on large scales it is sourced only by 
the nonlinear version of the entropy perturbation. We then expand these equations to second order 
in the perturbations, using a coordinate based formalism. Our results are relatively compact and 
elegant and enable one to identify the new effects coming from the non-canonical structure of the 
scalar fields Lagrangian. We also explain how to analyze, in our formalism, the interesting scenario 
of multifield Dirac-Born- Infold inflation. 
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1 Introduction 



Inflation offers a compelling mechanism to produce a spatially flat and approximately homogeneous universe. 
Moreover, it also provides a source for small primordial fluctuations that seed the observed large scale structures. 
Simple models of inflation are realized in terms of a single, slowly rolling scalar field characterized by canonical 
kinetic terms and a flat potential. In this case, the spectrum of fluctuations is characterized by adiabatic, almost 
scale invariant density perturbations, with nearly Gaussian distribution on super-horizon scales. 

While present-day observations are consistent with these models, there are good theoretical reasons to go 
beyond the simple hypotheses on which they are based. As an important example, string theory motivates 
frameworks in which many scalar fields, with non-canonical kinetic terms and not necessarily flat potential, 
play a role during inflation. Indeed, string theory predicts the existence of a large number of scalar fields, 
the moduli, which in the early universe may be sufficiently light to interact with the inflaton field, playing an 
active role during inflation. Moreover, the most studied stringy inflationary models are based on the dynamics 
of D-branes moving in higher dimensional spaces: in these set-ups, the inflaton field is governed by a Dirac- 
Born-Infcld (DBI) action characterized by non-canonical kinetic terms (see [1] for recent reviews). In some 
regimes, the non-canonical kinetic terms allow one to obtain inflation also with steep scalar potentials, as in 
the DBI- inflationary models of [2]- [5]. A more general analysis of the spectrum of primordial fluctuations for 
multi-field models with non-canonical kinetic terms, then, may allow one to understand how future observations 
will be able to probe, or exclude, inflationary models motivated or inspired by string theory. 

Single-field inflationary models with non-canonical kinetic terms, of a class belonging to models of ^-inflation 
[6] , have received much attention in the past years, for the possibility to generate spectra of fluctuations with non- 
Gaussian features, at a level that can be probed by future observations (see [7] for review and a comprehensive 
analysis of this case). More recently, models of inflation with more than one scalar field with non-canonical 
kinetic terms have been also considered [8]-[16]. In this case, non-standard correlations between adiabatic 
and entropy modes can occur, and this affects the amplitude of non-Gaussianities in these models [13]-[15]. 
These correlations can even persist in the primordial radiation dominated era, as first pointed out in [17]. The 
analysis of the examples studied so far suggests that the dynamics of fluctuations, in string inspired multi-field 
inflationary models, have new and rich features, depending on the form of kinetic terms for the inflaton field(s). 

In this paper, motivated by the previous arguments, we analyze cosmological fluctuations of a system of 
two scalar fields with non-canonical kinetic terms at a fully nonlinear level. We adopt the covariant formalism 
developed in [18]-[21], along the lines of earlier works by Ellis and Bruni [22] and Hawking [23] (see also [24]). 
This formalism is particularly suitable to analyze fluctuations of the system we consider, since it leads to a clear 
and natural decomposition of fluctuations in adiabatic and entropy components. This separation has been first 
studied, at the linear level, in [25] for the case of canonical kinetic terms (see [26] for a recent review). It is 
particularly convenient for analyzing the generation and conversion of adiabatic and isocurvature components 
of fluctuations produced during inflation. It also plays an important role in the previously mentioned recent 
papers, that analyze fluctuations for multi-field models with non-canonical kinetic terms. 

Following [21], we generalize this decomposition into adiabatic and isocurvature components to a fully 
nonlinear set-up. In order to do so, we first define adiabatic and isocurvature covectors, representing the 
nonlinear generalization of linear fluctuations, in the case of general field space metric. Then, we obtain a set of 
exact evolution equations for these quantities, in our framework characterised by a non-standard Lagrangian. In 
addition, by defining a suitable nonlinear generalization of the curvature perturbation, we derive the evolution 
equation for this quantity and we determine how it is sourced by the nonlinear entropy covector. This result 
generalizes already known linear equations that play an important role for analyzing how entropy fluctuations 
are converted into adiabatic modes on superhorizon scales. 

We then show how our results can be reexpressed in a more familiar coordinate based approach, and how 
our equations can be expanded to first and second order. Proceeding in this way, we re-obtain at first order 
the linear evolution equations determined in [11], while we find new results when pursuing the expansion up 
to second order. In particular, working in a large scale limit, we determine how the curvature perturbation is 
sourced by the entropy modes at second order. Our final equations are relatively compact, and allow one to 
clearly appreciate the impact of non-canonical kinetic terms on the evolution of fluctuations. 

Our discussion closely follows [21], to which we refer the reader for further details and references. We show 
how to extend the methods of this paper to the case of non-canonical kinetic terms for the scalar fields, and 
for general field space metric. In doing so, we obtain nonlinear equations for the fluctuations, that can be 
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expressed in a physically transparent way, and we show the usefulness of the covariant approach for analyzing 
cosmological models inspired by string theory. 



2 Covariant formalism for scalars with non-canonical kinetic terms 

We consider an arbitrary unit timelike vector u a — dx a jdr (u a u a = — 1), defining a congruence of cosmological 
observers. The spatial projection tensor orthogonal to the four-velocity u a is provided by 

h ab = g ab + u a u b , (h a b h b c = h a c , h b u b = 0). (1) 

To describe the time evolution, the covariant definition of the time derivative will be the Lie derivative with 
respect to u a , defined for a generic covector Y a by (see e.g. [27]) 

Y a = C u Y a = u b V b Y a + Y b V a u b , (2) 

and will be denoted by a dot. For scalar quantities, one simply has 

/ = u b V b f. (3) 

To describe perturbations in the covariant approach, we consider the projected covariant derivative orthog- 
onal to the four- velocity u a , denoted by D a . For a generic tensor, the definition is 

/u;,: : /<,/%: .../'/ ...v,//;;-. (4) 

In particular, when focussing on a scalar quantity /, this reduces to 

D a f = h b V b f = d a f + u a f. (5) 

We can also decompose 

V fc u a = cr ab + uj ab + -Qh ab - a a u b , (6) 

with the (symmetric) shear tensor a ab and the (antisymmetric) vorticity tensor to ab ; the volume expansion, 9, 
is defined by 

e = v a u a , (7) 

while 

a a = u b V b u a (8) 

is the acceleration vector. 

Let us now consider N scalar fields minimally coupled to gravity with general Lagrangian density (see [11] 
for a study of the linear perturbations in this type of model): 



Here X is 

X 

or, using Eq. (5), 



C = P{X,4> 1 ). (9) 



iGug^Va^V^ , (10) 



X = \ Gu^ 1 ^ - l -G u D a ^D a 4> J (11) 

and G[j = G[j(4> K ) is a metric in 'field space' that can be used to raise and lower field indices, denoted by 
capital letters. One can consider more general Lagrangians, like the ones studied in [13, 15]. 

However, the choice of Lagrangian density Eq. (9) is suitable for describing most of ^-inflationary models 
considered in the literature. Especially, it is important to point out that our formalism can be applied to 
two-field DBI inflationary models of the type studied in [12, 13]. One of the present author (S. RP) has shown 
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in these works that the multi-field DBI Lagrangian can not be written in the form P(X, (p 1 ). It can however be 
written in the form P{X , (j) 1 ) where X and X only differs by terms in spatial gradients. Therefore, all the results 
derived in this paper in the large scale limit, where spatial gradients can be neglected, are readily applicable to 
multi-field DBI models simply by considering the specific Lagrangian 



P(X, tf) = -y^jy (7l - 2fW)X lj Vtf) 



(12) 



where / is the so-called warp factor and V is a potential term (see [13] for more details on multi-field DBI 
models). 

The energy momentum tensor derived from Eq. (9) reads 

T ab - PxGuWa^Vb^ + g ab P . (13) 

Given an arbitrary unit timelike vector field u a , it is always possible to decompose the total energy momentum 
tensor as 

Tab = {p+ P)u a u b + q a U b + U a qb + 9ab P + ^ab, (14) 

where p, p, q a and ir a b are respectively the energy density, pressure, momentum and anisotropic stress tensor 
measured in the frame defined by u a . Starting from the expression for the energy-momentum tensor (13) one 
finds (we write <fti = Gu 4> J ) 

p = T ab u a u h = P x tffa - P , (15) 

P = \h ac T ab h b c = ^PxGuDatfi'D"^ + P, (16) 

q a = -u b T hc h c a = -PxfaDacp 1 , (17) 



ir ab = KT cd h d b -ph ab = Px (GuD^Dbtf -^-GuD^D^y 



(18) 



The evolution equations for the scalar fields are obtained from the variation of the action with respect to 
the fields themselves. In our case one gets 

Px GuV a V a <p J + Px T IJK (V a (f> J ) (V a <p K ) + (V a P x ) G u Va(f> J + Pj = (19) 

where Tijk = Gil^jk = \ {Gij,k + Gik.j — Gjk.i) is the Christoffel symbol associated to the metric G/j. 
The previous equation reduces to the usual one, in the case one chooses P = X — V{4> 1 ), where V is a potential. 

Notice that using the equalities, following from (5) and (6), 

D a D a = V a V V + ft + 6 - a b Drf 1 , (20) 

V a P x V a J - D a P x D a cb J - Px J , (21) 
V V V a <p K = D a (f> J D a (f) K -j> J j) K , (22) 

equation (19) becomes 

¥ + Fjk (<P J <P K ~ D a <j> J D a <j> K ) + {& + ^ tf-j^ G IJ Pj - D a D a ^ - a b D b ^ - D b tf = . 

(23) 

To simplify the notation, it is useful to define the spacetime derivative of field space vectors in curved coordinates 

VaA 1 = VaA 1 + T I JK V a <l> J A K , (24) 
with which we define a time derivative in field space 

V U A' = u a V a £ (25) 
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By choosing 



and a spatially projected derivative in field space 

n = K d K ...hf... v d Tlt- . (26) 

Notice that V a acts as an ordinary time derivative on field space scalars (i.e. quantities without field space 
indices) and T> a Gu = 0. Using these definitions, it is possible to rewrite (23) in a more condensed form as 

Vj 1 + Iq+^L) ^ - -^-G IJ Pj - V ±a (D> 7 ) - a b D h cp' - Drf 1 = 0. (27) 

3 Single field 

We start our discussion examining the simple situation of a single field. This provides the opportunity to 
introduce some quantities that will play an important role in the following discussion. 

Let us denote by <j> the single scalar field we are considering and set the trivial field space metric Gn = 1. 
Then the Klein-Gordon equation (27) reads 

D * Da <t> - a " D ^ - ^f'^ D b <f> = 0. (28) 
P.x dcp P.x 

Uc ° m ~ (29J 

one finds D a (j) = 0. Then (28) becomes formally identical to the well known homogeneous version, although it 
remains fully inhomogcncous and nonlinear. 

One can proceed further, and derive an exact and covariant equation that mimics the equation of motion 
governing the linear perturbations of a scalar field in a perturbed FLRW (Friedmann-Lemaitre-Robertson- 
Walker) spacetime. The idea is to consider the evolution of the space-time gradient of the scalar field, i.e., the 
covector 

0a = V a 0. (30) 

Indeed this quantity can be decomposed into a spatial gradient and a longitudinal component, 

(p a = -lj)U a + D a (f). (31) 

resembling, within our full nonlinear setting, the usual decomposition of a quantity in a homogeneous part, plus 
a perturbation. Notice that when making the particular choice (29) the spatial gradient disappears in the above 
expression, and indeed the evolution equation becomes in this case formally identical to the homogeneous one. 

We then construct a second-order (in time) evolution equation for <p a . We recall that a dot stands for the Lie 
derivative along u a , as defined in Eq. (2). We can derive the evolution equation for <j) a by taking the spacetime 
gradient of Eq. (28) and noting that, for any scalar <p, the Lie derivative with respect to u a and the spacetime 
gradient (but not the spatial gradient) commute, i.e., 

V a 0=(V a 0)\ (32) 

Acting with the space-time gradient V a on equation (28), and using the previous definitions, we obtain 



la+ I 9 + 5f ) ^ + < ^ Va 



= -<pV a e + V a (D b D b (b) + V a (a b D b ^)+V a 





■ i dp~ 


,x 


' 1 dP~ 




p^ 


~Px~ #. 



D b P 



x 



,x 



(33) 



This equation is similar to the analogous perturbation equation at linear order, but it additionally incorporates 
the fully nonlinear dynamics of the scalar field perturbation. 
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3.1 Integrated expansion perturbation on comoving slices 

One can define a covariant generalization of the comoving curvature perturbation using appropriate combinations 
of spatially projected gradients. For a scalar field, a natural choice is the covariant integrated expansion 
perturbation on comoving hypcrsurfaces 1Z a , defined as 

TZ a = -D a a + %D a <j>, (34) 

where a is the integrated volume expansion along u a , 

a ee i J drO (6 = 3d). (35) 

Since <d/3 corresponds to the local Hubble parameter, one sees that the quantity a can be interpreted as the 
number of e- folds measured along the world- line of a cosmological observer with four- velocity u a . 

When the four- velocity u a is chosen to be comoving with the scalar field, as defined in Eq. (29), then the 
last term in the definition (34) drops out. 

Note a useful property of TZ a - one can replace in its definition (34) the spatial gradients D a by partial or 
covariant derivatives, 

K a = -\7 a a+^ a . (36) 

9 

The same property applies to the nonlinear generalization of the comoving Sasaki-Mukhanov variable for a 
scalar field, which can be defined as 

Q a = D a <j>-tD a a=tn a . (37) 
a a 

3.2 Integrated expansion perturbation on uniform energy density slices 

Following [18, 19], it is also possible to generalize the curvature perturbation on uniform energy density hyper- 
surfaces. The key role is then played by the covector ( a defined as 

( a = D a a - -D a p. (38) 
P 

For a perfect fluid, the quantity ( a satisfies a simple first-order evolution equation 

e 2 

Ca = jrT a , (39) 

where 

T a = D a p-^D a p (40) 
P 

is the nonlinear nonadiabatic pressure perturbation (p is the pressure measured in the frame defined by u a , 
defined in Eq. (16)). For a barotropic fluid, r a = and £ a is conserved on all scales. The relation (39) for 
( a can be seen as a generalization of the familiar conservation law for £, the linear curvature perturbation on 
uniform energy hypersurfaces [28]. In the next sections, we will analyze the form of T a in various physically 
interesting situations. 

For a single scalar field, the comoving and uniform density integrated expansion perturbations ( a and lZ a 
satisfy 

(a+n a = --(D aP -^D a <f ) ) , (41) 

p V 4> J 

which simply follows from their respective definitions. The right hand side can be interpreted as the "shift" 
between hypersurfaces of constant p and hypersurfaces of constant cf) and the term inside the parenthesis 
represents the nonlinear generalization of the so-called comoving energy density perturbation of a single scalar 
field. 
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By choosing u a = u° om as defined in Eq. (29), the energy- momentum tensor of a single scalar field can be 
written in the perfect fluid form, i.e., with vanishing q a and 7r a j, and the energy density p and pressure p given 

by 

p = P X <j> 2 -P, p = P, (42) 

as can be checked by specializing Eqs. (15-18) to a single field and setting D a <f) = 0. 

After substitution in the definition (40) of the nonadiabatic pressure covector T a , we have 



P 



x 



1 



1 



Pa, + P 



D„ 



where we are using D a (f> = (implying also 2X — <p 2 ). We also used 

2 = E£ = P-X 
' ' P,x P x + xx ' 

Then, the evolution equation of ( a for a single scalar field is given by 



Ca 



D a <, 



3P; 



i+4 



P. 6 + P. 



(43) 



(44) 



(45) 



where we have used p = — P^x4> 2 (that can be obtain from the equation of motion for 4>) to get this expression. 

We will show in Sec. 5.2, as a particular case of a more general discussion, that the quantity on the right 
hand side of equation (45) vanishes on large scales, so that ( a is conserved in this limit. 



4 The two-field case 

4.1 Definition of adiabatic and entropy covectors 

In the two-field case, it is convenient to introduce a particular basis in the field space in which various field 
dependent quantities are decomposed into adiabatic and entropy components. In the linear theory, this decom- 
position was first introduced in [25] for two fields. For the multi-field case, it is discussed in [29] in the linear 
theory (see also [30]) and in [31] in the nonlinear context. 

In our case, the corresponding basis consists, in the two-dimensional field space, of a unit vector e£ defined 
in the direction of the velocity of the two fields, and thus tangent to the trajectory in field space, and of a unit 
vector el defined along the direction orthogonal to it (with respect to the field space metric) , namely 

ci = — , Gije# = l, Gueie J a =0, (46) 
o~ 

with 

a=^G u ^^ . (47) 

Notice that generically, the quantity a is not the derivative along u a of a scalar field a ; it is merely a notation. 
An important consequence of the above definitions is the identity 

S I j=c I rT e a j + c I s e sJ . (48) 

From e a ieJ a = 1, one deduces that is proportional to ef. It is then convenient to define 9 by 

V u ei = 9ei, V u ei = -0ei. (49) 

Again, it is simply a short hand notation; 9 is not the derivative along u a of an angle 9, although such an angle 
can be defined if the field space metric is trivial [21]. 

Making use of the basis (46), one can then introduce two linear combinations of the scalar field gradients 
and thus define two covectors, respectively denoted by a a and s a , as 

a a = e ff/ V a 7 , (50) 
s a = e s/ V a 7 . (51) 
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We will call these two covectors the adiabatic and entropy covectors, respectively, by analogy with the similar 
definitions in the linear context [11]. Whereas the entropy covector s a is orthogonal to the four-velocity u a , i.e., 
u a s a — 0, this is not the case for a a which contains a 'longitudinal' component: u a a a = a. It turns out to be 
useful to introduce the spatially projected version of (50-51), 

^ = e^iDacj) 1 =(T a + GU a , = c s/ D a 7 = s a . (52) 



4.1.1 Adiabatic Klein-Gordon equation 

We start our analysis of the evolution equations considering the adiabatic combination of the Klein-Gordon 
equations, i.e., the contraction of (27) by e„i. By noting that we can write 

a = eairf) 1 a = e„/Xi u f , (53) 

and by defining 3 

Pa = ciPj , (54) 

one obtains 



a + ( 9 + ^ ] a - -±-P, a - c aI V ±a (Z?V) - - <r£ - 0. (55) 

, X J , X , X 



The fourth term can be rewritten as 

e aI V ±a (£>>') = V ±a (c aI D a ^) - (V ±a c aI )D a ^ = D a a£ - (c sI V ±a ci)s a , (56) 
where we have used the definition of a^, Eq. (52), and the identity (48). Inserting 

Vj 1 = u»V b (V a /) + (V a u») (V^ 1 ) (57) 
into the calculation of the Lie derivative of s a with respect to u a , one gets 

esiVaei = -eaiVael = ^r(s a + 9<r a ) , (58) 
a 

from which one obtains 

e<Tl V ±a (£>°0 J ) = D a ai - i(« a + 6<j^)s a . (59) 
The adiabatic combination of the Klein-Gordon equations can thus be written as 

* + ( e + pf) & ~~Px~ = vv - _y W' (60) 



where we have defined 



and used the property 



Yw = v(i„ + e^)s a ^rf< > (6i) 



J}\H«X=VX, (62) 
which is valid for any covector orthogonal to u a (and thus also for s a ). Note that our coincides with the 
one of [21] in the case of canonical kinetic terms, in which situation the last term in (61) vanishes. 
Second derivatives of the scalar fields are hidden in the term P t x ■ Explicitely, this gives 



J + Q & + ^Z. = VV- - Y {s) + Z^L^^ + s a s a y , (63) 

c s Px Px *P,x 

where we have introduced the important, spacetime dependent, quantity 



P lip • (64) 

P,x + & Pxx 

Notice that c s is a fully nonlinear quantity; this implies for example that one cannot replace a 1 by 2X in its 
definition, since these quantities differ by terms in spatial gradients (see (11)). At linear order however, c s 
coincides with the now well known speed of sound of fc-inflationary models [6, 11]. 



3 Similar notations will be employed in the following, such as P Xa = eiPxi, P,xs = e£P,xj • 
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4.1.2 Entropic Klein-Gordon equation 

Let us now consider the entropic combination of the Klein-Gordon equations, i.e., the contraction of (23) with 
e s j. By using 

e sI V u ^ = 6&, (65) 
and by defining the entropic gradient of the potential, 

P,, = eiPj, (66) 

one finds 

Px 6-6 - P s - P x c sI V ±a (D a ( j> 1 ) - P x a a s a - D b P x s b = 0. (67) 
One can rewrite the last three terms by using the identity 

e sI V ±a (£>V) = D a s a + \{s a + 0a^ a , (68) 

derived by inserting (57) into the calculation of the Lie derivative of a a with respect to u a . Applying the 
property (62) to the covector s a , one finally gets 

&6-^- = V - s " + Y (°) > ( 69 ) 

with 

Y (a) ^-(^ + 0a^ a + ^f^ Sa . (70) 

0~ f y x 

To conclude, we have shown how to replace the Klein-Gordon equations for the fields by two new equations 
describing the evolution along the adiabatic and entropy directions respectively. Our covariant equations (63) 
and (69) look very similar to the homogeneous equations: on the other hand they capture the fully nonlinear 
dynamics of the scalar fields. Whereas these equations generalize the background evolution equations, we will 
proceed in the next subsection to derive fully nonlinear and exact equations for the covectors, which mimic and 
generalize the linearized equations for the adiabatic and entropy components. 

4.2 Evolution of the covectors 

We now derive evolution equations for the covectors a a and s a . More precisely, our purpose is to find two 
evolution equations, which are second order in time (with respect to the Lie derivative along u a ) and which 
mimic the equations obtained in the linear theory for the perturbations Sa and 5s (see [11]). 

4.2.1 Evolution equation for the adiabatic covector 

Starting from the definition of a a in (50), one finds that its time derivative, i.e., the Lie derivative with respect 
to u a , is given by 

&a = CalValj) 1 + 9.S a = V a & + 6s a , (71) 

where the last equality is obtained by using (j) 1 = ae^. A further time derivative yields 

<r a = V a <7 + 6s a + 6s a . (72) 
The next step consists in using (60) to eliminate a in the above expression. This gives 

»a - V a y (s) , (73) 



p.. 



oP 
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where we have used the relation 

V a P CT = P aa a a + P as s a + ^(s a + 8a a ) + P Xa V a X , (74) 

a 

and introduced the notation 

P aa = eieiVjVjP, P ss = eleiVjVjP, P sr7 = ^VjVjP . (75) 

In the previous expressions, T>i denotes the covariant derivative associated with G/j (thus having VjVjP = 
Pu-TFjPk). 

The evolution equation for a a can then be decomposed into a longitudinal part, obtained by contracting 
(73) with u a , and an orthogonal part obtained by contraction with /i ab . By using the relation 

Pa = P,**o + 0P S + Pxa X, (76) 

it is not difficult to see that the longitudinal part yields in fact the time derivative of (63). What is more 
interesting is the orthogonal or spatial part. Using several equations collected in the appendix, one can expand 
the various quantities inside the previous expression in terms of the perturbations and s a . After a long but 
straightforward calculation, the evolution equation for the adiabatic covector can be written in a form close its 
linear counterpart [11]: 



P 



(^r+le + iM^) +(c2-iRr 



Px 

+&c 2 s D a Q c 2 s D a (V c a c ^) 



a 

Z. 



-p. 



6P< 



crPxa)' 



P,XaZ( s ) 



H-H( 2 — 



Pxa 



M 



-c 2 s D a Y (s) + — (1 - c 2 s )D a U + -2- [P tXa - P XX Z {S) + e- 3a (e 3a aP xx y] D a U 



2a y ~ ~ s '-"- 2P X 
In the previous equation we introduced various quantities: 



(77) 



n = 



1 



Y, 



(s) 



aP 



,x 



[(1 + c 2 s )P s - c 2 s b 2 P. x& \ , 



(78) 
(79) 

(80) 



and 



M = 



2 P,sZ( s ) C 2 s Z( a) 



29 



,xx 



n 



■X 



p2 
r ,X 



,x 




2P X & 



P Xs - 



XX 



.X 



Moreover, we used the integrated volume expansion a introduced in Eq. (35). 



(81) 



Let us make a few comparisions with the case of a standard Lagrangian: the friction term is not only 
given by the volume expansion 9 but has a non vanishing contribution coming from the nonlinear speed of 
sound c s , which may vary in time (c s ^ 0) and may be different from unity (c 2 ^ l). The last two terms in 
(77), which exhibit a rich dependence on the form of the kinetic term, also vanish for a standard Langrangian 
P = X — V^ 1 ). Moreover, in that case, S reduces to 28, which encodes all the multi-field effects at linear order, 
as shown in [25]. The important parameter 5 plays an equivalent role for Lagrangians of the form P(X,<p I ) 
studied in this paper, as shown in [11]: if S = 0, the situation is effectively single-field at linear order. 

4.2.2 Evolution equation for the entropy covector 

Let us now consider the evolution equation for s a . From Eq. (58) the time derivative of s a is given by 

s a = o-e s /X> a e^ - 6a a . (82) 
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(83) 



Taking another time derivative, one finds 

s a = -Off a - Off a + °(s a + Off a ) + <xV J - a 2 R^ LJ e sI e J a e*V a <f> L , 
a 

where we have used (58) and (49) and the equality 

V a (VuA 1 ) - V u (VaA 1 ) = R I KLJ \7 a cf, L <j ) J A K (84) 

applied to the vector e£ and where R I KL] is the Riemann tensor associated to the metric G/j 4 . We can now 
use the entropic equation (69) to get rid of V a 0. Furthermore, using the relation 



V P S = P. sa ffa + PssSa -^(s a + 6<J a ) + P Xs V a X, 

a 



we obtain 



,x 



Sa--[ff- T r L )s a -[-^ + o 2 ) Sa -v a (v c s c ) 



p 



.X 



— — V a X + a R KLJ e sI e rT e a V a (f> 
P,x 



= -28a a P s + 



P 2 , 



x 



6,.. -Per, 



■X 



ffa + Va^(cr) • 



As for the adiabatic equation, the longitudinal part of this equation, upon using the relation 

P s = P„ s & - 0P. a + P Xs X, 
yields the time derivative of Eq. (69). The orthogonal part, instead, yields 



1 



P 



.x 



P 



X 



+ Z - a A R IKL jeteteie« \s a -D a (V c s c ) - -g± D a X 



.x 



= -26{a a t-P, s ^- 

1 V 



P P 

ff P,x P,x 



(ffaf+ D a Y [a) , 



(85) 



(86) 



(87) 



(88) 



where we have used the property that the covectors s a and s a are purely spatial, i.e., that (s a ) = s a and 
(s a ) 1 ' = s a . Using (69) as well as various expressions collected in the appendix, we finally obtain 



2P 



.x 



,Xs 



(ffa) 1 - + 
P,XxP,s 



U OC' t 

a I 5 _Zi 



Sa ~ D a (V C S C ) 



a \ c 



& 



K) ± + D a Y { 



,x 



with 



2P,XsPs 
p2 

,X 



(90) 

which mimics its linear counterpart in [11]. Let us remark that the last term in (89), absent in the linear theory, 
also vanishes for a standard Lagrangian. 

Starting from the fully nonlinear Klein-Gordon equations, we thus obtained a system of two coupled equations 
(77) and (89), controlling the evolution of our nonlinear adiabatic and entropy components. Although they 
appear quite involved, we will learn that they simplify considerably in physically interesting situations, in 
particular when considering linear or large scale limits. Moreover, as we will see later, immediately deduced 
from these equations are the already known evolution equations for the linear adiabatic and entropy components. 
Furthermore, since our equations are exact, they can be used to go beyond the linear order, up to second or 
higher orders in the expansion. 



4 We thus have R I KL] 



UL1 KJ KL LM KJ JM KL ' 
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4.3 Generalized covariant perturbations 

In this subsection, we will be first interested in the covariant generalization of the comoving energy density and 
curvature perturbations in the context of a two-field system. We will then consider the generalization of the 
curvature perturbation on uniform energy density hypersurfaces. 



4.3.1 Comoving energy density covector 

Let us introduce the covector 



D a p 



(91) 



which can be interpreted as a covariant generalization of the comoving energy density perturbation. In order 
to obtain the explicit expression of e a in terms of a a and s a , let us rewrite the components (15-18) of the 
energy-momentum tensor in the form 



P = P,x6- 2 - P, 
Px 

v = -fn + p, 

-Px&o-i, 



la 



= p,x ( n ab - -h ab W. 



where we have defined 

U ab = a^a^+ s a s b , 

while the quantity II has already been introduced in (79). 
Using (92), one finds 



(92) 

(93) 
(94) 

(95) 
(96) 



D a p 



( /> Y A <^ - P x ) [aL > a & - + 2P. X &D a & - P a a^ - P s S a + G 2 {Pxa^ + PxsSa) 



P = (Pxxv 2 - Px) I era - - + 2P x d-a - P> + P x * & 3 , 



(97) 



from which one obtains 



(98) 



This expression can be employed to rewrite Eq. (89) for the entropy covector s a in a useful different form. 
As noticed above, the longitudinal projection of Eq. (86) yields the time derivative of (69), which reads 



0_P^ 
d-Px 



Pr 



P Xs X P X P 



X *,s 



P 



X 



\xo 



aP 1 



(99) 



This relation, together with Eqs. (60) and (98), enables us to reexpress (89) as 



So + 



6 + 



Px Z (s) 
Px & 



Sa + 



ci a \ a 



s a - D a (X7 c s c ) 



< + D a Y (a) 



(100) 



As in the linear theory [11], this gives us an alternative expression for the evolution equation for s a , in which 
the comoving energy density perturbation appears explicitely on the right hand side. This expression will be 
useful in Sec. 5.2 when discussing the large-scale evolution of s a - 
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4.3.2 Comoving curvature covector 

Together with the comoving energy density, the comoving curvature perturbation can be generalized. For the 
general case of several scalar fields, this is done by defining the comoving integrated expansion perturbation 
(recall the expression for q a in (94)), 

Oi 

K a = -D a a- —^q a . (101) 

The definition of the Sasaki-Mukhanov covector given for a single scalar field in Eq. (37) can then be extended 
to the case of several fields, by defining for each field 

Qi = Datf - ^D a a . (102) 
Thus, the comoving covector lZ a can also be written as 

Ka = JJ^J- M a - (103) 

{4>J<t> ) 

In the iwo-field case, the definition (101) reduces, using (94), to 

TZ a = -D a a+-a^. (104) 
a 

Furthermore, one can generalize the Sasaki-Mukhanov variable to the adiabatic covector by defining 

Qa = C aI Qi =<?a- -Da®- (105) 

a 

Following [21], instead of giving the evolution equation of 1Z a , we will use a fluid description by considering the 
covariant generalization of the uniform density curvature perturbation, i.e., the integrated expansion perturba- 
tion on uniform density hypersurfaces £ a , defined in Eq. (38). In the two-field case, these two quantities are 
related by 

Oi 

(a+K a = --e a . (106) 

p 

In contrast with the case of a single scalar field, which can always be described as a perfect fluid, the total 
energy-momentum for two or more scalar fields corresponds in general to that of a dissipative fluid: the nonlinear 
formalism developed in [20] will thus be useful in this case. 

The adiabatic Klein-Gordon equation (63) can be rewritten as a continuity equation for the total energy 
density (92) and pressure (93), which reads 

p + Q(p + p)=V, (107) 

with the dissipative term 

v = p x (&z (s) + ien + ±nj . (ios) 

In [20] it was shown that the evolution equation for ( a for a dissipative fluid, which generalizes (39), is given 

by 

Ca = |t (T„ + S ) , (109) 

where the second source term on the right hand side, due to the dissipative nature of the fluid, is defined in 
terms of V as 

S «--q \D a V--D aP \ + —lD a e--D aP \ . (110) 

Let us return to the covector T a . Substituting the expression (93) for p into its definition (40), as well as 
using 

P - P> + Px (&a - ^tlj , (111) 

D a P = + P sSa + P.x (Wa - - \ D -^) > ( 112 ) 
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one obtains 



T a = -ef- [(1 + c 2 s )P,„ - c 2 s a 2 P XrT ] + P x *Zs a - P tX c 2 (dJ\ - % a pj + 1 (d u (P x H) - ^^-D a ^ 

(113) 

where we have introduced e a defined in Eq. (91). The above equation expresses the nonlinear nonadiabatic 
pressure perturbation in the two-field case, which sources Eq. (109). It will be useful in the next section, where 
we will consider both the linear and super-Hubble approximations of our evolution equations. 

The evolution of £ a , governed by Eq. (109), is thus sourced by a rather complicated term obtained by 
summing T a and S a in Eq. (110). We will learn in the next sections that this term simplifies considerably in 
physically interesting cases, when either the linear or the super-Hubble approximations are taken. 



5 Approximate equations 

In this section we study the evolution equations of the two-field system under two types of approximations: 
the linear limit and the limit in which we neglect higher orders in spatial gradients. In an expanding FLRW 
universe, the latter corresponds to the large scale limit. In the rest of the paper we will use the symbol ~ to 
denote an equality at the linear level, and the symbol w to denote an equality valid only on large scales. 



5.1 Homogeneous and linearized equations 

In many cosmological applications, since our Universe appears to be close to a FLRW universe on large scales, 
it is sufficient to restrict oneself to the linearized version of the evolution equations. We first consider this 
linearization procedure directly at the level of the covariant equations, as done in [22] . 
In a strictly FLRW universe, all the spatial gradients vanish and therefore 

a£=0 , s a = . (FLRW) (114) 

Consequently, the scalar quantities Y/ s \ and Yt a \, defined respectively in (61) and (70), vanish and the evolution 
equations for a and 9, respectively (60) and (69), reduce to 

ct+ (zH+^) a-^f- = 0, (FLRW) (115) 

= (FLRW) (116) 

where we have introduced the Hubble parameter H — 6/3. Not surprisingly, the above equations exactly 
correspond to the homogeneous equations given in [11]. Furthermore, in a FLRW universe, all the terms in 
Eq. (77) for erf and Eq. (100) for s a vanish. 

At linearized order, we treat the covectors erf and s a , vanishing at zcroth order, as first-order quantities. 
Similarly, their derivatives s a , s a , (<r ) and (<j a ) 1 ' are first-order quantities. Therefore the linearized version 
of the evolution equations are simply obtained by keeping only the homogeneous terms in the coefficients 
multiplying the spatial projection of a a , s a and their derivatives. The linearized versions of (77) and (100) are 
thus (recall the definition of S in (80)) 




(*.) X + 



3H + 



P 



x 



Px V c 2 s 
+&c 2 s D a Q - c 2 D a (D c a^) 



Px 



-P** - Px8 2 + e- 3a (e 3a aP, Xa y 



^ (Hc 2 y p c 



Hcl 



,X<7 



H 
H 



and 



s a + { SH + 5^ ] s a +(^ s + =^) Sa -D a (D c s c ) 



,x 



P,x° 



(117) 



(118) 
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Note that the terms involving Y( s ), F( CT ) and II have disappeared, since these scalars are quadratic in first-order 
quantities. We have also replaced V c <7^- by D c a£, as well as V c s c by D c s c , since their difference is quadratic 
in first-order quantities, according to (62). Indeed, the acceleration vector a b , which vanishes at zeroth order, 
is considered as a first-order quantity. 

One can also linearize the evolution equation for £ a . As discussed above, the terms containing and II 
can be neglected. The dissipative term V thus reduces to 

T>~P,x&D a a£, (119) 

while the expression for e a becomes 

e a - - Sot - *3 S „] ■ (120) 

On the other hand, the expression for T a becomes 



T a * -e a j [(1 + c 2 s )P„ - cl& 2 P X a] + PxcrEs a , (121) 



from which it is immediate to obtain the expression for £ a . 



5.2 Expansion in spatial gradients 

Apart from the linearization procedure, there is another approximation in the cosmological context that is used 
to describe the Universe on very large scales. This approximation is based on an expansion in spatial gradients, 
which are small for scales larger than the local Hubble radius [32, 33, 34]. Recently for instance, it has been 
adopted in [35] for studying nonlinear perturbations by means of the ADM formalism, in a system describing 
a single scalar field with arbitrary kinetic terms, up to second order in the gradient expansion. In the present 
article, we restrict ourselves to the leading order: in this perspective one sees, from their definition (52), that 
and s a are first-order quantities with respect to spatial gradients because they are linear combinations of spatial 
gradients. The scalars and Y/^ however are second-order with respect to spatial gradients since they are 
quadratic in and s a (or their time derivatives). Hence, the right hand side of Eq. (60) and of Eq. (69) can 
be neglected on large scales, so that these two equations become, in the large-scale limit, 

(122) 



(123) 




9- 



0. 



Although they look very similar to the homogeneous equations (115) and (116), these equations are fully 
inhomogeneous and encode the evolution of nonlinearities on large scales. This limit illustrates the separate 
universe picture [36, 37] where the inhomogeneous universe can be described, on large scales, as juxtaposed 
Friedmann homogeneous universes. 

If, so far, the order in spatial gradients seems to coincide with the perturbative classification of the previous 
subsection, it differs however for the term 



(124) 



which is first order perturbatively but second order in spatial gradients, at least for the first term on the right 
hand side, since a£ is already first order in spatial gradients. For the second term, it has been shown in [21] 
that u a can be chosen so that a c is at least first order in spatial gradients. We will explicitly verify it in the 
next section by working in a coordinate system. 

With these prescriptions, the evolution equation of <r^ and s a obtained at lowest order in spatial gradients 
become 



(*«)" 
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,x 



P<T 



PxO 



- 3a (e 3Q crP 



Xa 



at + <Jc 2 s D a Q 



,xcr 



(125) 
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and 

-p-y e « ( 126 ) 

where we have dropped the terms containing II in (100), which are two orders higher than s a in spatial gradients. 

We now expand the evolution equation for £ a , Eq. (109), by neglecting higher-order spatial gradients in the 
two terms on the right hand side of this equation. The nonadiabatic pressure perturbation becomes 

T a « -e a j [(1 + c 2 s )P a - c 2 s a 2 P Xa ] + P, x &Sa a , (127) 

while the dissipative nonadiabatic pressure perturbation E can be completely dropped, since the dissipative 
term V is at least second order in the spatial gradients and thus S a is third order in the spatial gradients. 
Equation (109) therefore becomes, on large scales, 



Ca « [(1 + c l) p ,° - & 2 P,x*] e a - ^s a . (128) 



Note that the lowest order limit in spatial gradients of the evolution equations of s a and £ a , respectively 
Eqs. (125-126) and (128), are similar to their linear counterparts except for the terms D a {D c ajr) and D a (D c s c ), 
which are third order in spatial gradients and therefore negligible in the spatial gradient expansion. This is 
because in these equations the terms that are higher than linear order in the perturbative expansion turn out 
to be also higher than first order in the spatial gradient expansion. 

We now concentrate our attention on the comoving energy density perturbation, e a , defined in Eq. (91). 
Until now we have only made use of the Klein-Gordon equations of the scalar fields. However, in order to study 
the behavior of the comoving energy density, we will now make use of the Einstein equations, in particular of the 
so-called constraint equations. The projection of the Einstein equations along u a yields the energy constraint in 
a covariant form, 

u a G ab u b = 8ttG P . (129) 

Assuming that the vector field u a is hypersurfacc orthogonal, it is possible to use the Gauss-Codazzi equations 
and the decomposition 

D b u a = a ah + -Qh ab , (130) 

which is the spatially projected version of (6) (with u ab = since u a is here hypersurface orthogonal), in order 
to rewrite the energy constraint as 

I (^) R + le 2 - a ab a ab ^j = SirGp, (131) 

where is the intrinsic Ricci scalar of the space- like hypersurfaces orthogonal to u a . 
The mixed projection of Einstein's equations yields the covariant momentum constraint 

u h G bc h c a = 8irGq a , (132) 

which can be rewritten, via Gauss-Codazzi relations and Eq. (130), as 

D b a a b - ^D a e = 8irGq a . (133) 

By combining the energy and momentum constraints, one obtains the nonlinear covariant version of the general- 
ized Poisson equation, which in the linear theory relates the comoving energy density to the Bardeen's potential 
defined from the curvature perturbation. Here one finds 

X -D a {^R - a bc a bc ) + &D b o- b a = 8nG e a , (134) 

where we have introduced on the right hand side the quantity 

i a = D a p-eq a = D a p + P x e&ai. (135) 
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This quantity can be seen as an alternative generalization of the comoving energy density perturbation, since, in 
the linear limit, it is equivalent to e a defined in (91). In the fully nonlinear case, the two quantities in principle 
are different: using (107) one gets 

la - e„ = * - lp,x®Il) °a- (136) 



3 

This difference becomes however negligible on large scales, 

e a ~ e a . (137) 

Now, the left hand side of Eq. (134) contains the projected gradient of the Ricci scalar, Dj^R. From its 
definition in terms of derivatives of the metric, it can be shown that, for a perturbed FLRW universe, this term 
is of third order in the spatial gradients. Equation (134) can thus be used to show that e a in Eqs. (106), (126) 
and (128) can be neglected on large scales, if the shear can also be neglected in this limit. Indeed, on large 
scales, the shear rapidly decreases in an expanding perturbed FLRW universe. Thus, in this limit the comoving 
and uniform density integrated expansion perturbations ( a and lZ a coincide (up to a sign), 

(a+TZa^O. (138) 

Furthermore, one can rewrite Eqs. (126) and (128) as a closed coupled system of equations, describing the 
large-scale nonlinear evolution of adiabatic and entropy perturbations, 

s„ + ( B f ^ j .s„ : ( H 2 S + ^\s a ^Q (139) 




and 



Ca « -|^«a. (WO) 



Remarkably, these equations for the covectors £ a and s a look very similar to their linear counterparts, studied 
in [11], although they extend them to describe the full nonlinear evolution for the perturbations. Despite their 
compact appearence, they exhibit a quite involved dependence on the form of the Lagrangian. For instance, 
the mass term ji s , which reduces to V ss — 2 for standard Lagrangian, also depends in our case on the entropic 
projection of the mixed derivative P,xi, a characteristic shared by the important parameter 5. 

As explained in Section 2, we stress again that the equations derived in this subsection in the large scale 
limit are also valid for two-field DBI inflation with the adequate Lagrangian Eq. (12). In the next sections, 
always following [21], we will show how to translate our results, obtained within the nonlinear formalism, in 
the more familiar coordinate-based approach, obtaining quite straightforwardly the evolution equations for the 
perturbations up to second order. 



6 Linear perturbations 

In this and next section, we relate the covariant approach to the more familiar coordinate based formalism. We 
analyze the linear perturbations in the present section, and we then consider an expansion up to second order 
in the next one. 

Let us thus introduce generic coordinates x^ — {t, x 1 } to describe an almost-FLRW spacetime. Here a prime 
will denote a partial derivative with respect to the cosmic time t, i.e. ' = d/dt, since the dot has been used till 
now to denote the Lie derivative with respect to u a . 

The background spacetime is a FLRW spacetime, endowed with the metric 

ds 2 = g^dx^dx" = -dt 2 + a{t) 2 lij dx i dx j . (141) 
At linear order, the spacetime geometry is described by the perturbed metric 

ds 2 = (g^ + 5g^)dx»dx v , (142) 
where the components of the metric perturbations can be written as 

Sgoo = -2 A Sg 0i =aB l , 8g l3 =a 2 B ir (143) 
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We decompose, as usual, Bi and Hij in the forms 

B l = W l B + BY, (144) 
= -2ifa ij +2Vrf j E + 2$ (i EY ) +2E%, (145) 

where BY and EY are transverse, i.e., ViB Vt = = WiE Vt , and Efj is transverse and traceless, i.e., V ' iE T% i = 

and -f^Efj = 0. Here Vj denotes the three-dimensional covariant derivative with respect to the homogeneous 
spatial metric 7^ (which is also used to lower or raise the spatial indices). The matter fields are similarly 
decomposed into a background and a perturbed part, 

<f> I (t,x i ) = 4?(t)+6<f> I (t,x i ). (146) 

We now need to specify the components of the unit vector u a , which defines the time derivation in our 
covariant approach. At zcroth order, it is natural to take it orthogonal to the homogeneous slices. At first 
order, [21] chose, for simplicity, such that Uj = 0. This implies that, up to first order, the components of 
are given by 

u" = {1 - A, -B l /a}, (147) 
and those of the acceleration vector are given by 

a" = {0,VM/« 2 } . (148) 

This confirms that can be considered as first order in spatial gradients, in agreement with the assumption 
of the previous section. 

Since the formalism relies on many covectors, it is useful to first consider a generic covector Y a and work 
out the components of its time derivative Y a . To make the explicit calculation, it is convenient to replace in the 
definition of the Lie derivative (2) the covariant derivatives by partial derivatives and write 

Y a = u b d b Y a + Y b d a u b . (149) 

At zeroth order, the components of Y a are simply 

^ M = {y„'. 0,0,0}, (150) 

assuming that the spatial components Yj vanish so as to respect the symmetries of the geometry. At first order, 
we get from (149) and (147) 

8(Y ) = 6Y£ - (YqA)' , S(Yi) = 5Y! - Y d t A. (151) 

Let us now consider, for the special case of two scalar fields, the adiabatic and entropic covectors <r a and s a , 
which we have introduced earlier. Note that with the choice of four-velocity (147) ay- = <7j. The background 
equations of motion can be deduced immediately from Eqs. (63) and (69) and read 

%+3Ha' + ^a> 2 ~^k = 0, (152) 

C-s P ,X P.X 



a'9' - 

Px 

with 



(153) 



a> = yjGufi'p'. (154) 



Note that we used the symbol 6' to represent the background value of 9. Again, as there is no angle 9, this is 
merely a notational convenience. 

From its definition, Eq. (50), one finds that the spatial components of <r a at linear order can be expressed as 

Sat — eaidid^ 1 — dida . (155) 
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We use the notation 



Sa = eaidiScfr 1 , 



(156) 



in agreement with previous studies in the literature. By using Eq. (151) with (154-155) one finds that the 
spatial components of the first and second time derivatives are given by 



5(&i) = d t (Sa' - a' A) , S{a l ) = d l (Sa" - a' A' - 2a" A) . 
The same procedure for s a gives 

Ssi — di5s, Ss = esidiScj) 1 , 



(157) 
(158) 

which also coincides with the notation of [11]. Since s a , in contrast to a a , has no longitudinal component, s = 
and the spatial components of s a and s a are simply 



S(si) = Bids' , 6(§i) = diSs". 



(159) 



Plugging the explicit components (157) and (159) into the linearized equations for a a and s a , given by (117) 
and (118) respectively, one easily obtains the linearized equations for Sa and Ss. These read, respectively, 



Sa" + 



3H + 





c 2 r 
Sa' + JX- 




r.x L 



-P,aa - PxO 2 + e~ 3a (e 3 "a'P x *)' 



Sa 



A' + I 2"- + 3H (1 + c 2 ) + jX- ) ) A + 3c 2 V' - c 2 s V 2 (E> B/a) 



%V 2 Sa 
or 



= ESs' + 



E'-E 2^- 



C S P. 



and 



Ss" +[3H+^\Ss'+[fi 2 s + ^)Ss- -^V 2 8s 



Ss 



n2 

r,2 i " 



P,x* 



where we have used 



6 = 3H, 8Q = -3HA — 3tp' + X7 (E 1 — B /a) 



(160) 



(161) 



(162) 



In the latter equation we have introduced the first-order comoving energy density perturbation Se, defined 



by 



Sti = diSe, Se = Sp — —Sa, 

a' 



(163) 



which follows from the definition (91) of e a . Using 

P 



P.xa' 2 -P. 



P.x a> 



Sp = ~^-(Sa' -a'A-ZSs) + (Pxaa' 2 -P a )8a, 



one sees that Se can be expressed as 

Se = 



Px a> 



Sa' — a' A -Sa — ^Ss 



Moreover, linearizing the spatial components of the energy constraint (131) yields 

3H (i)' + HA) - -^V 2 [tp + H(a 2 E' - aB)} = -4ttGS P , 

while the momentum constraint (133) gives, since Sqi = —P.x a'diSa, 

ip' + HA = AttGP x aSa. 



(164) 



(165) 



(166) 



(167) 
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Combining the two above constraints yields the relativistic Poisson-like equation 

1 



rV 2 [V> + H(a 2 E' - aB)} = AnG5e, 



(168) 



which can also be directly obtained by linearizing the spatial components of Eq. (134). This equation shows 
that the comoving energy density perturbation 5e is second order in the spatial gradients, and thus negligible 
on large scales in Eq. (161). 

The quantity 5a is not gauge-invariant, in contrast with 5s. This is why it is useful to consider the gauge 
invariant Sasaki-Mukhanov variable Qsm, defined as [38] 



Q SM = 5a + —tp. 



(169) 



Note that the above traditional definition does not follow exactly from the definition of Q a given earlier in 
Eq. (105). Indeed, from Q a , one can extract a scalar quantity Q defined as 



Qt = d t Q, Q = 5a- -5a, (170) 
where 5a can be written in terms of metric perturbations by making use of Eqs. (35) and (162) (see [19]), 

5a = -ip+^J V 2 (E' - B/a)dt. (171) 



Thus the scalar variable Q coincides with Qsm only in the large-scale limit. 
In the flat gauge 5 , defined by $ = 0, 5a coincides with Qsm, 



5a = Q 



SM- 



(172) 



In this gauge, it is possible to use the momentum constraint equations (167) to derive the metric perturbation 
A as a function of Qsm, 

i = -^7Q S M, (173) 



and one can write the Poisson equation (168) as 
V 2 (E'-B/a) = --^-„_ 



Ha' 



Q'sm + 



JT~~a 7 



Qsm - ^5s 



(174) 



where we have used the expression (165) specialized to the flat gauge. 

By replacing Eqs. (173) and (174) into the evolution equation of 5a, one finds the evolution equation of Qsm 
[11], ' 



Qsm+|3H + ^-(^# 



with 



r 2 V 2 \ / (He 2 ) P 

Qsm = (^5s) - s 



a* 



He 2 



a'P 



x 



(°'in)' 

a'/H 



3H 



c 
~P 



.X 



(a'/H)' \ (a'/H)' 
a'/H ) a'/H ' 



5s, (175) 



(176) 



When one considers only large scales, the expression (174) reduces to 

H' 



Q'sm + ( -jj ~ -rj ) Qsm -^Ss^O, 



(177) 



which means that there exists a first integral for the quantity Qsm and that the second-order equation of motion 
(175) is not necessary in this limit. In fact, one can easily check that the large-scale limit of (175) is an automatic 
consequence of the first integral (177). 



5 In the following, a hat indicates that the corresponding quantity is evaluated in the flat gauge. 
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Let us now consider the evolution equation for ( a with two scalar fields. The spatial components of ( a , at 
linear order, are given by [19] 

Sd = diC, ( = Sa- ^5p. (178) 

From Eq. (171) the scalar variable ( is thus related to the Bardeen gauge invariant variable (b, defined as 
[39, 40] 

Cb = -1> - ^Sp, (179) 
P 

in a similar way to how Q is related to Qsm, and on large scales these two quantities coincide. 
According to Eq. (151), the spatial components of ( a are 

S(Ci)=diC'- (180) 

On large scales, one then finds the relation 

C^-^ESs. (181) 
a 

Note also that from Eq. (106) there is a simple relation between £ and 72, 

( + ll = ~^-5e, (182) 
P 

which shows that £ and —72 coincide on large scales. 

In conclusion, all the equations in this section, derived directly from the covariant formalism, exactly repro- 
duce the linear results of [11]. In particular, according to Eqs. (161) and (175), whereas entropic perturbations 
are sensitive to the usual Hubble horizon (k w aH), adiabatic perturbations are amplified at sound horizon 
crossing, c s k rts aH. As shown in [12, 13], these equations are not directly applicable to two-field DBI inflation, 
whose Lagrangian is not of the form studied in this paper. Despite that, in that situation, it turns out that 
Eq. (175) for the adiabatic perturbation is still valid and that Eq. (161) for the entropy perturbation is mini- 

mally modified, simply replacing the gradient term -^V 2 Ss by ^§-V 2 (5s, for the adequate Lagrangian Eq. (12), 
with the consequence that both types of perturbations are amplified at sound horizon crossing. 

In the next section we turn to the second-order perturbations, and show the power of the method we are 
using to derive new results. 



7 Second order perturbations 

In this section we consider second-order perturbations (for a recent review, see [41]). We will thus decompose 
any scalar quantity Y as 

Y(t, x l ) = Y(t) + (t, x l ) + SY^ (t, x l ), (183) 

where Y(t) is the background part and SY^ and SY^ are respectively the first and second-order contributions 
(note that we do not follow here the convention of including a numerical factor 1/2 in front of the second-order 
contribution). In our subsequent equations, to simplify the notation, we will often omit the index W for the 
first-order quantities, unless it is required for clarity reasons. 

Our main purpose will be to expand our equations governing a a , s a and ( a at second order in the perturba- 
tions. We start with a brief discussion on the gauge invariance of the curvature perturbation. As shown in [19], 
the second-order expression of the spatial components Q can be written, after some manipulations, in the form 



Sp 

with 



Cf=^C (2) + |5 4 C (1) ' (184) 



C (2) = 5a (2) _ H Sp( 2) Sp 
P P 



(185) 



and being given in Eq. (178). In [21] it has been shown that is gauge invariant on large scales: their 
proof holds identical in our case and we then refer the reader to this paper for more details. We now pass 
to discuss the evolution equations for adiabatic and entropy fields up to second order, and analyze how they 
source the curvature perturbation on large scales. The reader should keep in mind that the following discussion, 
restricted to large scales, is applicable to two-field DBI models on using the adequate Lagrangian Eq. (12). 
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7.1 Adiabatic and entropy fields 



Here we derive the evolution equations for the adiabatic and entropy field perturbations at second order. We 
will restrict ourselves to large scales, a situation in which the equations become more tractable, starting from 
the equations expanded in spatial gradients discussed in Sec. 5.2. For convenience, we have collected in the 
appendix various background and first-order expressions that will be used in the rest of this section. 

Second-order evolution for the perturbations when several scalars are involved has been first studied in detail 
by Malik in [42] in the large-scale limit and with canonical kinetic terms (see also [43, 44, 45]), while it has been 
analyzed in [46, 47] using the separate universe approach. The second-order decomposition into adiabatic and 
entropy components however, has only been derived in [21], in the case of standard kinetic terms (see [48] for 
a recent application and [49, 50] for related studies). This decomposition can easily be derived since the fully 
nonlinear adiabatic and entropy components have already been identified. We will show how to extend this 
decomposition into adiabatic and entropy components to the case of non-trivial metric in field space, and use 
it to derive the corresponding evolution equations for these modes. 

In order to pursue the decomposition, one starts by expanding the definition of a a in Eq. (50) at second 
order. After some straightforward manipulations, one can write the spatial components of <7j as 



5a? ] = a 



SWtf + t-etfjjK (e*e{5a 2 



e«e{6s 2 



2e^eidaSs) 



+ ^7 (5s' + 9' 6a) d l 5s. (186) 



To deal with the term Ss'dids, which cannot be written as a total gradient, it is convenient to introduce the 
spatial vector 



Vi = ^{Ssd.Ss' - Ss'diSs) = -5s' 8,5s + ^d^SsSs'), 

which vanishes when 5s' and 5s have the same spatial dependence, i.e., 5s — f(t)5s. 

By expanding also the definition of s a in Eq. (51), one finds, for and <7j, respectively, 



(187) 



5af ] = d t Sa (2) + ^5o-di5s - IjV 



a' 
5a. 



5s™ = dM 2) + °^5s\ 



(188) 
(189) 



5a^ = ew^V + \eit JIK (e^So* + ef e^s 2 + 2ef e[SaSs) + ^j5s5s', 
5 8 ™ = e a7 <j(V + \e J s t JIK (ee> 2 + ef e{5s 2 + 2e^e I s 5a5s) - % (ss' + ^5a 



with 

(190) 
(191) 

The form of the right hand side of Eq. (189) has been chosen by analogy with the form (184). Since s a vanishes 
at zeroth order, arguments similar to those used in [21] to prove gauge invariance for C (2) ensure that 5s ( ' , 
defined in Eq. (191), is also gauge invariant on large scales. The form of 5s 1 - 2 - 1 is, then, forced by our covariant 
definition. Let us remark that the non-flat nature of the field space metric manifests itsef in the appearance of the 
terms with Christoffel symbols in 5a^ and 5s^ 2 \ Otherwise, these are formally identical to the corresponding 
definitions in [21]. 

(2) 

Note that 5s\ contains the first-order adiabatic perturbation. This is due to the fact that the adiabatic 
and entropy components are defined locally: whereas the first-order components are defined with respect to a 
background basis in field space, which is only time dependent, the second-order components will be sensitive to 
the first-order fluctuations of the field space basis, which can be expressed in terms of the first-order adiabatic 
and entropy components. Then, the adiabatic component a a does not vanish at zeroth order and 5a^ is not a 
gauge invariant variable. It is then useful to consider our generalization of the Sasaki-Mukhanov variable, Q a , 
defined in Eq. (105). Its spatial components can be expanded at second order in the perturbations, similarly to 
what we have done with ( a , 

Qf ) = d,Q^ + ^Q' + ^Qd t 5s - IjVt, (192) 



a' 
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where is defined as 



q(2) = 5a m °' ( 2) Sa 
H H 



From this expression it is natural to define 



(193) 



(194) 



as the local part of the scalar gauge invariant second-order Sasaki-Mukhanov variable. Restricted to a single 
scalar field, this definition coincides with the one given in [42]. Note that we cannot write Eq. (192) in the 
same form as Eq. (184) because the last two terms on the right hand side cannot be written as a total spatial 
gradient. 

The second-order (in time) evolution of a a is given by Eq. (125). However, on large scales we do not need to 
compute a second-order differential equation because the adiabatic evolution is governed by a first integral, as 
in the linear case. This first integral is obtained directly from the constraint equations and it is not necessary 
to expand (125) at second order in the perturbations. 

In order to compute this first integral we need the second-order energy and momentum constraints, which 
can be derived by expanding Eqs. (131) and (133) and by using 



~ ^ HA 2 + 3A ^> _ Q^> _ 3HA (2) _ 3^(2)' _ 
On large scales, one can write the second-order energy constraint equation as [21] 



3H 



HA m + ^(2)' + 2A(HA + V/) 

2 H 



-4ttGV 2) , 



(195) 



(196) 



where 5p^ can be calculated using (note that the non-flat nature of the field space metric simply manifests 
itsef in the presence of the Ricmann tensor) 



5^ a 2 = a' 2 (AA 2 - 2A^) - AAa' (5a' - 1 5 s) + (5a' - fSs) + (Ss' + 1 5 a) ' 
- 2W , e'5s^ + 2a'5a {2) ' - 26" ' 5a Us' + ®-5a j - a' [^5s5s'^j 



+ a l2 eieiefe L s R KIJL 5 S 2 . 



Explicitcly, this gives 



5 P 



(2) 



P 



X 



(a'SaW - a' 2 A^ - a'ESs^) + (a' 2 P, Xa - P a ) 5a^ + A p 



where A p is a quadratic function of A, 5a^ x \ 5s^\ and their first derivatives, 



A 



= i§5a'{5a' 29' 5s) + *L 



± ((1 + c 2 )P a - c 2 s a' 2 Pxa) 5s - (P s - a' 2 P. Xs ) 



5a 



+ 



P.aa - v' 2 Px** + j, (i% - 1 ' ' Px s) 5a 2 + (a' 2 P Xaa - P,„) 5s5a 



a' 2 P.Xss - Pes + % ( nl 



+ 6" + a' 2 eieiefe L s R KIJL 



5s 2 



+ \ (Pxxx^ + 3P tXX a' 2 ) (5a' - d' 5s - Aa' 

(P.xxvv' 2 + Pxv) Sa + (P XXs a' 2 + P Xs ) 



+ a 



5s A 



(5a' - e'5s - Aa') . 



(197) 



(198) 



(199) 
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The second-order momentum constraint equation reads [21] 



HA^I + ^ (2) ' + 2W - ^HA 2 - A(HA + i/,') 



(2) 

where the second-order momentum 5q^ is given, from Eq. (94), by 



-inGSq, 



(2) 



Sq^ = -ft 



.X 



P.x o 

2c 2 a 



7 5a 2 + P iX 0'5a5s 



—5ediSa + PxVi. 
a' 



(200) 



(201) 



(2) 

As already noticed in [31] for large scales, 5q\ cannot be written as a total gradient when several scalar fields 
are present. After neglecting 5e on large scales in the above equation, this is manifest because of the presence of 
Vi. This implies that, in principle, if Vi does not vanish, one cannot dehne at second-order a comoving gauge, 
i.e., such that 5qf^ = and 5qf^ = 0, in contrast with the linear theory or the single-field case. 

However, it is instructive to derive the evolution equation for Vi on large scales by using the linear evolution 
equation for 5s, Eq. (161), neglecting the gradient term and 5e at first order. One finds a simple expression 
that generalizes the corresponding result of [21] to the case of non-canonical kinetic terms 



v;+ ( 3 #+5f J Vi = o. 



(202) 



This implies that, in an expanding universe, P^xVi will decay like a~ 3 and rapidly become negligible even if it is 
nonzero initially. Consequently, in an expanding universe one can in practice ignore Vi on large scales and thus 
define, in an approximate sense, a comoving gauge at second order, which coincides with 5a^ = = Sa^ 2 \ In 
this approximate comoving gauge, the momentum (201) can be written as a total gradient. In the rest of the 
paper, in order to remain as general as possible, we will keep the term Vi. 

Similarly to the first-order case, it is possible to combine the energy and momentum constraint equations 
and derive the relativistic Poisson-like equation analogous to Eq. (168), which corresponds to the expansion, at 
second order and on large scales, of Eq. (134). By expanding Eq. (135) and using (137), one has 



5e[ 2) w d l 8 P {2) - ZH5q (2) - SSSq, w 0, 



(203) 



where the last approximate equality is a consequence of Eqs. (162), (167), (196) and (200) and confirms our 
conclusion of Sec. 5.2 in a covariant context, namely that we can neglect e a on large scales. The second-order 
spatial components of e a defined in Eq. (91), can be decomposed as 



a' 



5a. 



with Se^ defined by 



5e^ = 5p 



(2) 



-8a 



(2) 



5a 



5eW 



IE 

2\a' 



5a 



9'5s 



(204) 



(205) 



It is only when Vi is negligible that the quantity 5e^ can be interpreted as the comoving energy density at 
second order. Otherwise, as discussed before, the comoving gauge cannot be defined. 

Using the decomposition (204) and the fact that 5e^ is negligible on large scales, Eq. (203) can be written 

as 

d 2 8e^ w ZHPxd^i- (206) 

When Vi is negligible, and only then, one finds that, like at first order, the second-order comoving energy density 
is negligible on large scales, 

<5e (2) w 0. (207) 
Having discussed the general properties of the second-order constraint equations, we can now derive the 

(2) 

evolution equation of the gauge- invariant adiabatic component Qg^- Similarly to what we have shown in the 
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(2) 

previous section for the first-order variables, the simplest way to derive an equation satisfied by is to work 

,(2) 
*SM> 

(2) „(2) 



in the flat gauge -0^ = = ^ 2 \ where 5a^ reduces to 



6a 



ISM- 



In this gauge, Eqs. (196) and (200) reduce, respectively, to 



<3H 2 (i( 2 ) - 2i 2 ) w -AnGSp 



(2) 



Hdi { i(2) _ |i2 



a ™ (2) 



Using the first-order constraint equations Eqs. (173) and (177), these can be rewritten as 



i(2) 



1 



and 



3H 2 + H'/c 2 
H' 



2c 2 . 



- Z6s^) + i (P CT - ^Px.) Qsm + 6 (fr)" QIm - ^ 



A, 



i( 2 ) « - 



Ha' 



)(2) 
iSM 



-Qsm<5s + — 
r 2cr 



3- 



P 



-.d~ 2 d l V l 



(208) 

(209) 
(210) 

(211) 
(212) 



This last equation contains a nonlocal term because we have written the momentum constraint as a scalar 
equation while keeping the second-order vector Vi defined in (187). By combining these two relations to get rid 
of AW one obtains the following first integral for Qg^: 



o {2) ' + 



EL _ £1 \ n( 2 ) _ ^A S (2) 



Sa'Pxcl - — ( _ - 3- 



H' 



H 



<5sM 



- \ZHc 2 + 



H' 



-QsmSs - -d- 2 d*Vi 



(213) 



This equation is the second-order equivalent of (177). We have put together on the left hand side all the terms 
depending on purely second-order quantities, and on the right hand side all the terms which are quadratic in 
first-order quantities. As in the first-order case, the entropy perturbation sources the evolution of the adiabatic 
perturbation. The nonlocal term containing V. comes from the momentum constraint and is a new feature with 
respect to the first-order case (or the second-order case for a single scalar field). However, as we have discussed 
earlier, it becomes quickly negligible on large scales in an expanding universe, in which case the first integral 
(213) becomes a scalar local equation. 

We eventually derive a second-order equation for the entropy mode Ss^ 2 \ which will be second order in time. 
Since we are restricting ourselves to large scales, we simply expand the spatial components of Eq. (126) up to 
second order. This gives 



5(Si) 




(2) 



-.2 , " 



.,(2) 



+ S\ G 




+6 I M 2 + ^ 



(2) 



Pxa'~ 1 ' 



(214) 



where we have neglected the gradient of the comoving energy density at first order, 5ep which according to 
Eq. (168) is subdominant on large scales. 

To proceed, we need the spatial components of the first and second time derivatives of the covectors s a . 
By using Eq. (149) for s, at second order and ignoring the higher-order terms in the gradient expansion, one 
obtains 

«5(s,) (2) ^58^' -ABM, (215) 
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and, by applying once more (149), 

«5(s,) (2) « SaW" - A'd l 5s' - 2Ad l 5s" . (216) 

Recalling that 

5sf } = dM 2) + ^Ss' , 

we can then substitute the expressions (189) and (215-216) for the second-order entropy component and its 
derivatives, and use equation (161) and its time derivative. Defining 



-X '° ° & 2 P X 



M u = + & 2 R IKJL e^e L a + [(2 + c 2 )PjPj + c 2 & 4 P XI P XJ - c 2 s b 2 P XI Pj - c^PxjPj] 

(217) 



such that 

j2 



+ — = e 1 s e J s M IJ , (218) 
the equation (214) can be written as the spatial gradient of the following scalar equation 

5s {2) " + Uh + SaW + (ft + J) Ss^ « - 6 -,5s' 2 (zH^ |e^ s J M /7 ) W 



ele J s r 



<r'(S - 0"') M IJiX + ef %M ;J -2§M„ 



Ss 2 - -=^=7(5e (2) . (219) 



Let us make a few clarifying comments: when calculating the various derivatives of Mjj with respect to the 
fields and to X, b 2 is replaced by 2X and the basis vectors e£,e£ are considered as constant (they arc not 
functions of X and the fields); eventually, T> K M U = d K M u - t^jM LJ - t^jM IL . 

Note that this equation is closed, in the sense that only the entropy field perturbation appears: even when 
Se^ is not negligible, it can be written in terms of Ss and 5s' by using Eqs. (206) and (187). Thus, on large 
scales the entropy field evolves independently of the adiabatic components, as in the linear theory. 



7.2 Generalized uniform density and comoving curvature perturbations 

We now derive the large-scale evolution equation for £( 2 ), expanding to second order Eq. (128). In the fluid 
description, it was shown in [51, 19] that on large scales the evolution equation for can be written as 

C (2)'~ ^_ r ( 2 ) L_ nCi, (220) 

^ p + P p + P V ; 

where can be read from the second-order decomposition of the quantity r a for a fluid, defined in Eq. (40), 
i.e., 

5if> -a 4 r( 2 ) + ^a 4 r (1) '. (221) 

In the two-scalar field case considered here we must compare this expansion with the expression for T a given in 
the large-scale limit by Eq. (127). Expanding this equation to second order, using Eq. (189), one obtains 

Sr\ 2) « A (2) {^[(l + ^)F,-^]} + (^'S) (diSsW + yd^ 

+ d l SsS(P x o-E) (222) 
where we have used that Se^ is negligible on large scales. From this, defining £j as 

£/ = (1 + c 2 )Pj - c 2 s a 2 P XI (223) 

such that 

P x bZ = e I s Z I , (224) 
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one obtains the expression 



eie^j.j-a'ie - H)e& 



.Y 



2 



(225) 



From the previous expression it is immediate to rewrite the second-order evolution equation for (( 2 ) on large 
scales, Eq. (220), as 



C 



(2)/ 



H 



'x <?'■ 



1 



p 2(7 



a 



S)e^-2PxH 2 ) ^ 



6^ 
2 

(226) 

The previous expression shows that the second-order curvature perturbation is, as expected, sourced only by 
entropy modes on large scales. On the other hand, it depends in a much richer way on £/ and on the derivatives 
of P, that is on the form of the kinetic terms, than its linear counterpart in eq. (181). Moreover, reduces to 
— 2Vj for a standard Lagrangian, with no dependence on the kinetic term X. The appearence of £jx on the 
right hand side of (226), acting as a source for the large-scale evolution of the curvature perturbation, is thus 
solely due to the non-standard nature of the Lagrangian we are considering. It is important to emphasize that 
the nonlinear formalism we adopted, as well as the fact that we considered a very general Lagrangian, have 
allowed us to straightforwardly obtain a relatively elegant and compact expression. The same would be difficult 
to get working directly in a coordinate based approach, or considering a specific Lagrangian. 

It is also useful to express our results in terms of lZ a . The spatial components of lZ a can be decomposed as 



n f) = din v) + fo d . n v>-IL Vu 

a' <j' z 



with 



TZ^ ee -Sa^ + 5fcr (2) + S 4 



2 \ & I a' 



(227) 
(228) 



The last term in Eq. (227) comes from the fact that, like e a and in contrast to £ a , 7Z a is defined in terms of the 
spatial momentum which cannot be expressed in general as a pure gradient. When this term can be neglected, 
and only then, IZ^ coincides with the second-order comoving curvature perturbation defined in [52, 53]. 

It is easy to derive a first-order (in time) evolution equation for 7Z^ by noting that and 7Z^ are related 
on large scales. Indeed, expanding Eq. (138) to second order using (184) and (227), and neglecting terms 
proportional to Se^, one gets 

1 



C (2) +7e (2) 



3P 



,xcr 



12 



(229) 



When Se^ is negligible on large scales, like in an expanding universe where we can neglect Vi, and IZ^ 
coincide on large scales as in the single-field case [53]. However, this is not true in general in the multi-field case 
if Vi cannot be neglected. 

From this relation and the evolution equation of (( 2 \ Eq. (226), one can find a large-scale evolution equation 
forft< 2 \ 
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2a 



|£ 6sSs > + Ss W + (el^j o>0 ~ H)e&,* - 2P X S 2 ) ^ + (c 2 + ^)<J. 



H' 



c(2) 



(230) 

The second-order uniform adiabatic field perturbation TZ^ can be related on large scales to Q 1 ^, by combining 
Eqs. (194) and (228). One obtains 



TZ^ « ? 



4m - {Q'su - 6' 8s) Qsm - ( £ ) Q 



1 fH 



rr 



2H V a 



(231) 



which can be used, together with the linear first integral (177), to show that Eq. (230) is equivalent to the first 
integral (213). 
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8 Conclusions 



In this work, we adopted a covariant formalism to derive exact evolution equations for nonlinear perturbations, 
in a universe dominated by two scalar fields. These scalar fields are characterized by non-canonical kinetic terms 
and an arbitrary field space metric, a situation typically encountered in inflationary models inspired by string 
theory. Our exact equations can be expressed in a relatively compact way, and in physically interesting limits 
they closely resemble their linear counterpart. They acquire a quite simple and elegant form, due to the fact 
that we considered a very general Lagrangian to perform our calculations. 

Extending the methods of [21] to a general field space metric, we have decomposed our non-linear scalar 
perturbations into adiabatic and entropy modes, corresponding to a generalization of analogous definitions 
adopted in the linear theory. Then, we derived the corresponding evolution equations, that acquire several new 
contributions associated with the non-canonical kinetic terms for the scalar fields. We also obtained a nonlinear 
generalization of the curvature perturbation on uniform energy density hypersurfaces, showing that on large 
scales it is sourced only by the nonlinear entropy perturbation. 

We have used these nonlinear equations as a starting point to show how to extend and generalize some 
results previously obtained in the literature, in a relatively straightforward way. In particular, we went beyond 
previous works by computing explicitly the evolution of the second-order adiabatic and entropy components on 
large scales, in the case in which the scalar fields have non-canonical kinetic terms. The second-order adiabatic 
component is, in this limit, governed by a local first-order (in time) evolution equation, sourced by terms depend- 
ing on the second-order entropy perturbation as well as, quadratically, on the first-order entropy perturbation. 
Both first and second order entropy perturbations satisfy a second order (in time) evolution equation and the 
full system of equations, valid on large scales, is thus closed. The final system of equations enables one to 
clearly identify new effects due to the non-canonical structure of the scalar fields Lagrangian. A representative 
example is equation (226), that governs the second-order evolution equation for the curvature perturbation on 
large scales. This quantity is sourced by combinations of first and second-order entropy perturbations, with 
coefficients depending on the form of the kinetic terms. As discussed in Sec. 7.2, in some situations, cancella- 
tions can occur, reducing the size of some coefficients with respect to the others. In these cases, the curvature 
pertubation would be sensitive only to some type of terms, with important consequences for the evolution of 
second-order perturbations and thus for non-Gaussianitics. All the results derived on large scales in Sec. 7 are 
applicable to the interesting scenario of multifield DBI inflation of the type studied in [12, 13]. 

In conclusions, we have shown that our results, and more in general the covariant approach to the dynamics 
of fluctuations, can be important to analyze the evolution of perturbations in multi-field models of inflation 
inspired by high-energy physics or string theory. This interesting subject will offer, in the near future, unique 
opportunities to allow comparisons between predictions of high-energy physics motivated inflationary models, 
and observations of the cosmic microwave background radiation. 
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A Useful identities in a two-field system 
A.l Covariant identities 



Acknowledgments 



D a P, x 



D a X 



& a + GU a - 0S a = &a ~ 8s a , 

6- D a 6- - -^DaR, 

P,xx D a X + P Xa ui + P Xs s, 



(232) 
(233) 
(234) 
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D a P XX 
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+ P.Xa <J 



D a P X = (D a Pxx) 

+ P,xx 



P,XXX D a X + P^XXa a a + P.XXs s a ■ 
P.XXa D a X + P.Xcrcr <J~t + P.Xas S a , 
1 •' 



oa- -II 
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A. 2 Background identities 



{D a a)a--D a U 



H' = -47rGiW 2 . 



+ D a (P.Xa &) 

+ P.XXO- (7a - 0s a - 0S a 



(235) 
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(237) 



(238) 



(239) 



P[ a = a' (P itra + a"P X a) + 0'P S , (240) 

P; s = a' (P,„ + a"P Xs ) ~ O'P.a, (241) 

P[ as = a' (P, aas + o"P X <rs) + 0' (Pes ~ P.crcr) , (242) 

P; ss = a' (P CTSS + a"P X ss) - 26'Pvs, (243) 

P'x = o> (Px* + <t"Pxx) , (244) 

P'xx = ° (Pxxa + a'Pxxx) , (245) 

P [xa = °' (Pxaa + °" P,XXv) + 0'P,Xs, (246) 

P'xs = ° (Pxas + a'Pxxs) - e'Pxa, (247) 

(248) 

0' = (249) 

0" = ^ (Pre + d»P Xa ) J^rj (P„ + P,x^' 2 + d"Zf) , (250) 

P,x P.x<y' V c i ) 

A. 3 First-order identities 

Sei = 1 (5s' + 6' 5a) e{ - f ^4 (e«5a + e«6s) , (251) 

Sei = ~ (5s' + 6' So) ei T^ej (e* 5a + e«5s), (252) 

6{&) = 5a' - 6' 5s - a' A , (253) 

5P S = P sa 5a + P, ss 5s + P Xs (°'5{a) - ^5Uj - ^-(5s' + 6' 5a), (254) 

5P,x s = PxsaSa + PxssSs + P^xxs (&'5(o) - ^6Uj - ^{5s' + 6' 5a), (255) 

5P SS = P ssa 5a + P, sss 5s + Pxss (o'5{a) - ^5Uj - 2^ (5s' + 9' 5a) , (256) 

5P,x = PxaSa + P. Xs 5s + P X x (e'5(o) - ^5Uj , (257) 

5P X - (SPxY-APx- (258) 
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